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Abstract 


The  elanstic  scattering  of  eleartrons  by  atomic  hydrogen  including  the 
effects  of  exchange  and  polarization  is  calculated.  To  include  the  exchange 
effect  the  total  spin  and  space  function  of  the  incident  and  the  bound 
ekctrons  is  made  antisymmetric  with  respect  to  interchange  of  the  two  elec- 
trons. In  this  way  an  eqiiation  for  the  wave  fur:ction  of  the  incident  electron 
similar  to  the  Fock  equation  is  obtained.  To  include  the  effect  of  polarization 
of  the  hydrogen  atom  by  the  incident  electron  the  wave  function  for  the  bound 
electron  is  taken  in  such  a  way  that  the  electrostatic  potential  energy  between 
the  incident  electron  and  the  atom  with  this  wave  function  becomes  equal  to 
the  average  Coulomb  and  dipole  potential  energy  of  the  incident  electron, 
the  dipole  potential  being  determined  by  the  experimental  value  of  the  polar- 
izability  of  the  hydrogen  atom.  Since  the  exchange  potential  depends  on  the 
form  of  the  bound  electron  wave  function,  the  effect  of  polarization  on  the 
exchange  effect  is  also  taken  into  account.  The  equation  for  the  wave  function 
of  the  incident  electron  is  reduced  to  a  set  of  radial  integro-differential 
equations  by  the  method  of  partial  waves.  These  equations  are  integrated 
numerically  using  for  the  bound  electron  wave  function:  first,  the  ground 
state  hydrogenic  wave  function,  and  second,  the  wave  function  calculated 
according  to  the  above  scheme.  The  result  is  that  in  the  second  case  the 
calculated  curve  for  the  total  elastic  cross  section  reproduces  most  features 
of  the  experimental  curve  for  the  total  elastic  scattering  cross  section  of 
electrons  by  hydrogen  atoms  obtained  by  B.  Bederson. 

The  calciilation  in  the  first  case  has  been  done  by  P.M.  Morse  and  W.P. 
Allis.  It  seems,  however,  that  there  are  some  errors  in  their  formula. 
These  corrections  have  been  included  in  the  calculated  curves. 
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I.  Introduction 

mie  scattering  of  electrons  by  hydrogen  atoms  has  long  teen  vsed 
as  a  model  for  the  scattering  of  electrons  by  atomic  systems,  because 
it  is  amenable  to  some  sort  of  approximate  analysis  and  because  it 
exhibits  many  of  the  feattires  of  the  latter  problem.  In  order  that  the 
approximate  analysis  be  useful  in  the  more  complicated  systems,  some 
basis  must  be  fotmd  for  evaluating  the  accuracy  of  the  approximation 
techniques  used   in  this  elementary  problem.  To  that  end  an  extensive 
series  of  experiments,  involving  the  elastic  scattering  of  electrons  by 
hydrogen  atoms  in  the  range  between  zero  and  the  first  excitation  potential 
of  the  hydrogen  atom  ,  has  been  \mdertaken  by  Professor  Benjamin  Bederson 
and  his  collaborators  at  New  York  University'-  -■ ,  and  by  Dr.  Wade  L.  Fite 
at  General  Atomics  *-  -^ . 

Unfortunately,  the  resvilts  of  these  experiments  do  not  agree  with  one 
another,  the  principal  area   of  disagreement  being  the  fact  that  Prof. 

Bederson  finds  that  at  about  5  electron  volts  there  seems  to  be  a 

2 
resonance  in  the  cross  section  which  is  of  magnitude  lOOita  .  Fite,  at 

General  Atomics,  does  not  find  such  a  peak. 

Since  the  theoretical  investigations  to  date  have  found  no  peak,   one 

would  think  that  the  Fite  results  are  to  be  preferred  over  those  of 

Bederson.  The  three  most  important  types  of  calculations  which  have  been 

done  in  order  to  investigate  this  problem  have  been  a  variational 

[221 
calculation'-  -■ ,  and  a  calculation  using  a  Hartree  type  of  approach  including 

a  polarization  term"-  -' ,  and  one  involving  an  adiabatic  approximation'-  -• .   If 

the  Bederson  effect  is  correct,  and  a  true  resonance  exists,  it  may  not 

be  surprising  that  a  variational  calculation  would  not  uncover  it.  One 


wotild  expect  that  the  tirLal  function  used,   in  the  variational  principle 
would  have  to  resemble  the  correct  wave  function  to  a  considerable  degree 
in  order  for  the  resonance  to  appear  in  this  calciilation.  The  ordinary 
Hartree,  or  Hartree-Fock  method,  as  applied  to  this  problem,  does  not, 
as  will  be  indicated  below,  accurately  describe  the  polarization  of  the 
atom  by  the  incoming  electron.  In  the  calculations  that  have  been  done 
using  this  model,  one  can  correct  for  this  by  inserting  a  phenomenological 
polarization  term  into  the  Hamiltonian.  One  of  the  troubles  with  this 
method  of  attack  is  that  the  phenomenological  polarization  term  which 
one  inserts  mvist  be  cut  off  as  one  approaches  the  origin  of  coordinates, 
and  that  the  cross  section  is  very  sensitive  to  this  cut-off.  In  the 
third  method"-  -'  the  distortion  of  the  atom  is  accounted  for  in  i;he 
adiabatic  limit.  This  calculation  leads  to  an  extra  potential  term 
representing  the  effect  of  the  polarized  atom  on  the  incident  electron 
and  when  the  incident  electron  is  at  large  distance  from  the  atom  thJLs 
term  approaches  the  classical  value  for  the  potential  field  due  to  an 
induced  dipole.  The  validity  of  the  second  order  perturbation 
calculation  is  based  on  the  assirarptions  that  first,  the  incident  electron 
moves  slowly  eno\igh  30  that,  according  to  the  adiabatic  principle,  the 
excited  states  of  the  hydrogen  atom  ai*e  eigenfunctions  of  the  unperturbed 
hydrogen  atom,  and  second,  that  the  interaction  energy  between  the 
incident  and  the  bound  electron  is  small  in  comparison  to  the  absolute 
value  of  the  energy  of  the  ground  state  of  the  hydrogen  atom.  The 
validity  of  these  gissumptions,  especially  the  second  one,  may  be  questioned 
since  the  main  contribution  to  the  total  cross  section  is  from  the  region 
in  which  the  interaction  energy  between  the  two  electrons  is  not  small  in 

comparison  to  the  absolute  value  of  the  energy  of  the  ground  state  of  the 
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hydrogen  atom.  Thus,  the  absence  of  the  resonance, when  this  type  of 
calculation  is  used, is  not  necessarily  a  reason  for  discarding  the 
experimental  resvilts  of  Professor  Bederson. 

In  this  paper,  ve  shall  propose  a  model  for  the  scattering  of  electrons 
by  hydrogen  atoms,  which  includes  an  alternative  method  of  taking  into 
account  the  effects  of  polarization.   It  predicts  a  resonance  at 
approximately  the  same  energy  that  Bederson  found  and  the  cross  section 
is  of  the  same  order  of  magnitude. 

We  make  no  claim  that  our  method  of  including  polarization  effects 
is  superior  to  previoxas  attempts.  However,  when  one  considers  the  results 
obtained  in  this  paper  one  can  see  that  the  theoretical  models  are  still 
too  crude  to  enable  us   to  distinguish  between  the  two  sets  of  experimental 
data.  Consequently  our  conclusion  must  be  that  the  siibject  is  still  not 
exhaiosted.   It  is  to  be  emphasized  that  neither  oior  method  nor  any  of  the 
previous  methods  which  have  been  used  for  the  analysis  of  the  scattering 
of  electrons  by  hydrogen  atoms  can  make  any  estimate  of  what  the  errors 
involved  in  the  calculation  are  and  that  consequently  the  only  recourse 
to  the  correctness  of  the  method  is  comparison  with  the  experimental  results. 

One  of  the  first  attempts  to  analyze  this  problem  has  been  via 
a  Hartree-Fock  calcvilation,  done  by  Allis  and  Morse  "^  -*  .  The 
difficulty  in  choosing  a  product-type  wave  function  which  is  character- 
istic of  the  Hartree  or  the  Hartree-Fock  method,  is  that  no  self- 
consistent  proced\ire  is  possible  once  one  has  made  the  assumption 


* 

There  are  a  few  errors  in  this  paper.  In  the  second  term  of  the 

2  2 

right-hand  side  of  (5)  the  factor  k  -E  should  be  k  +E.  This  changes 

the  factor  k  -1  in  (8)  to  k  +1.  Fomrula  (9)  should  be  the  same  as  (8) 

with  the  right-hand  side  multiplied  by  -1. 
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that  a  product  wave  function  is  to  be  used.  The  reason  for  this  is  that 
in  order  to  satisfy  the  boundary  conditions  at  infinity  it  is  necessary 
that  one  of  the  terms  of  the  product  be  the  ground  state  hydrogenic 
wave  function.  Thus  one  loses  the  value  of  the  self-consistent  procedure 
in  that  one  of  the  terms  of  the  product  wave  function  is  now  fully 
determined.  However,  if  one  of  the  product  wave  functions  is  the  ground 
state  hydrogen  atom,  then  clearly  we  have  not  taken  into  account  the 
effect  of  the  distortion  of  the  atom  by  the  incoming  electron.  For  low 
energy  collisions,  this  is  a  serious  matter. 

We  have  taken  the  attitude  that  it  is  less  important  to  satisfy  the 
boundary  conditions  at  infinity  than  it  is  to  describe  accmrately  the 
effect  of  the  interaction  of  the  incoming  electron  with  the  hydrogen  atom. 
We  describe  the  effect  of  the  polarization  of  the  atom  by  the  incoming 

electron  by  saying  that  the  effective  charge  distribution  of  the  atomic 

-■yr 

system  in  a  static  description  is  given  by  e  '  ,  where  7  <  1  in  atomic 

units.  The  effect,  then,  of  the  polarization  is  to  spread  out  the 
charge  distribution  of  the  atomic  system.  7  is  a  parameter  which  is  to 
be  determined;  it  would  probably  be  different  for  different  angular 
distribution  of  the  incoming  electron.   It  may  even  be  a  function  of  the 
energy  of  the  incoming  electron,  although  we  have  not  used  such  a 
refined  model.  This  effective  charge  distribution  for  the  bound  electron 
is  now  substituted  in  a  Hartree,  or  Hartree-Fock,  type  of  calculation, 
and  the  res\iltant  integro-differential  equation  is  integrated  numerically 
on  the  lEM  lOk   machine. 

If  one  wished  one  could  satisfy  the  boundary  condition  at  infinity  by 
changing  the  effective  atomic  potential  when  r^  is  large.  We  do  this  by 


adding  to  the  potential  obtained  with  our  charge  distribution  another,  such 
that  the  sum  is  equivalent  at  this  distance  to  the  one  obtained  when  we 
use  the  ground  state  hydrogenic  wave  function.  Since  the  effective  potential 
for  r,  large  is  small  in  any  case,  it  is  not  expected  that  this  correction 
would  greatly  affect  our  results. 

In  the  following  section  we  discuss  the  general  technique  for  reducing 
the  three  body  problem  to  a  one  body  problem  and  the  way  in  which  one  can 
estimate  the  value  of  y   to  be  used.  In  Section  3  we  derive  the  equations 
of  motion  for  the  scattered  electron  on  the  basis  of  the  model  we  use. 
These  equations  are  integrated  in  Section  h.     The  results  and  conclusions 
are  presented  in  Sections  5  and  6. 

2.  General  considerations 

2.1  Reduction  of  the  three-body  problem 


Any  three-body  problem  can  be  reduced  to  a  two-body  problem  by  expressing 
the  two  independent  relative  position  vectors  in  terms  of  the  three 
position  vectors.  In  the  present  problem,  since  the  mass  of  the  proton 
is  large  coiapared  to  the  mass  of  the  two  electrons,  the  motion  of  the 
proton  can  be  disregarded.  It  is  then  convenient  to  take  the  position  of 
the  proton  as  the  origin  of  the  coordinate  system.  The  two  relative  position 
vectors  would  then  be  the  position  vectors  of  the  incident  and  the  bound 
electrons  with  respect  to  the  proton.  We  call  these  vectors  r^  and  r_, 
respectively. 

To  further  reduce  the  two-body  problem  to  a  one-body  .problem  we  assxane 
that  the  total  wave  function  of  the  incident  and  the  bound  electrons  is 
the  product  of  the  one -electron  wave  functions, 

+  ^  -A      +  -i  +  _^ 

(2.1)      ■^'{t^,x^    =  p-(r^)  ^^{r^   +  r{r^(}i{^x^   , 


where  ^"(r, ^r^)  is  the  total  vave  function  and  F  (r.)  and  0(rp)  are  the 
vave  fionctions  of  the  incident  and  the  bound  electrons,  respectively. 
The  upper  and  the  lower  signs  on  F  (r^ )  represent  the  two  distinct  states 
of  the  scattering  system  in  which  the  two  electrons  have  antiparallel  or 
parallel  spins,  respectively.  The  second  term  on  the  right-hand  side  of 
(2.1)  is  to  make  the  total  spin  and  space  wave  function  of  the  scattering 
system  antisymnetric. 

By  assuming  a  known  form  for  the  bound  electron  wave  function  0(rp), 
the  coordinate  r^   in  the  equation  for  the  two-body  problem  can  be 
eliminated  by  integrating  in  this  equation  over  coordinate  t^.     The  details 
of  this  integration  will  be  given  In  the  next  section.  In  this  w^  the 
three-body  problem  is  reduced  to  a  one-body  problem  and  the  resulting 
equation  which  depends  only  upon  the  coordinate  r  describes  the 
scattering  of  the  incident  electron. 

The  choice  of  0(r„)  is,  however,  important  because  0(rp)  determines 
the  effective  interaction  between  the  two  electrons;  in  turn  this 
interaction  determines  the  phase  shift  and  the  elastic  cross  section  of 
the  incident  electron. 

In  this  paper  0{r^)   is  chosen  so  that  it  closely  approximates  the  region 
of  interaction  of  the  two  electrons.  This  choice  is  based  on  the  following 
eirguments : 

From  a  classical  point  of  view,  when  the  incident  electron  approaches 
the  hydrogen  atom  it  separates  the  center  of  mass  of  the  negative  charge 
distribution  from  the  positive,  thereby  increasing  the  energy  of  the  atom. 
Since  this  separation  is  along  r^ ,  the  wave  function  of  the  botind  electron 
should  have  an  axis  of  symmetry  in  this  direction.  The  wave  theory  gives, 
however,  a  different  picture  of  the  charge  distribution  of  the  boiand 
electron.  This  can  be  understood  by  expanding  the  incident  electron  wave 
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(2.2)     F'ir^)-^Ii    ^/""i^^ji^^i^     '     ^^l  =  cos  0^. 


function  in  terms  of  the  zonal  harmonics, 

+  .      T   00    + 

■1  ^=0 

It  can  be  seen  from  this  equation  that  corresponding  to  each  angular 

momentum  quantum  nimiber  J?   there  is  a  definite  charge  distribution  of 

the  incident  electron.  In  particular,  when  J?  =  0,   the  charge  distribution 

of  the  incident  electron  is  spherically  symmetric  with  respect  to  the 

origin.  The  charge  distribution  of  the  incident  electron  determines, 

in  t\im,    the  charge  distribution  of  the  bound  electron.      In  this  way  we 

expect  that  for  .*  =  0  the  charge  distribution  of  the  bound  electron  is 

spherically  symmetric  with  respect  to  the  origin.     For  simplification  in 

the  wave  equation  of  the  incident  electron  we  have  assimied  that  this   charge 

distribution  is  spherically  symmetric  for  aU  values  of  ^  .     A13J.6  and 

Morse  made  precisely  the  same  assimptions  in  their  work.     The  simplest 

form,  for  such  a  wave  function  would  be  one  similar  to  the  wave  function 

of  the  grotind  state  of  the  hydrogen  atom*-  -■ , 

7^r 


(2.3)  0  (r)    =  Ne     °    ,      7     =—     , 

*  o  o  >oa 


N  being  the  nonnalization  constant,  and  a  ,  the  Bohr  radius.  To  account 
for  the  polarization  of  the  atom,  the  wave  function  of  the  perturbed 
hydrogen  atom  is  represented  by 

(2.1;)     0(r)  =  Ne"'^ 

with  y  <  y   .     In  what  follows  the  relationship  between  y   and  the  polarizability 
of  the  hydrogen  atom  through  a  phenomenological  arg\jment  is  shown. 

2.2  Connection  between  y  and  the  polarizability 

If  0(r)  is  to  specify  the  state  of  the  polarized  atom  it  should  produce 


an  electrostatic  potential  which  is  the  sum  of  the  Coiilomb  and  the  dipole 
potentials  of  such  an  atom.  As  far  as  the  relationship  between  7  and 
the  polarizability  is  concerned,  we  can  adopt  the  classical  point  of  view 
and  imagine  that  the  center  of  mass  of  the  negative  charge  is  displaced 
a  distance  d  from  its  unperturbed  position  in  the  presence  of  the  electric 
field  of  the  incident  electron,   d  is  connected  to  the  polarizability 
of  the  hydrogen  atom  through  the  relationship 

(2.5)  ed  =  p  =  a£      , 

where  p  is  the  induced  electric  dipole  of  the  hydrogen  atom,  a  its 

polarizability,  and  E  the  electric  field  of  the  incident  electron.  Since 

?  2 

E  =  -  e/rf^,  by  (2.5),  d  =  -  a/r^ .  When  the  incident  electron  is  at  large 

distances  from  the  atom  a  force,  due  to  the  induced  dipole  equal  to 

2pe/r^  =  -2a;e  /r^,  is  exerted  on  the  incident  electron  and  this  force  is 

derivable  from  a  dipole  potential  given  by 

k 

(2.6)  V^  =  ae/2rj^  . 

When  r  is-  small  the  main  electrostatic  interaction  between  the  incident 
electron  and  the  hydrogen  atom  is  the  Coulomb  potential  of  the  nuclexis 
given  by 

(2.7)  v^  =  e/r^  . 

We  can  derive  an  equation  which  reduces  to  (2.6)  and  (2.7)  for 
appropriate  values  of  r^  by  making  use   of  a  crude  model  of  the  hydrogen 
atom  in  an  electric  field.  Assume  that  when  the  hydrogen  atom  is  in  an 
electric  field,  the  electronic  charge  distribution  remains  spherically 
symmetric  but  that  its  center  of  mass  moves  a  distance  d  to  the  right  as 
is  shown  in  Figiore  1.  The  electric  field  is  due  to  an  electron. 


Figure  1 


-y   r" 
Since  the  hydrogen  atom  wave  function  is  given  by  0(r"  )  =  N  e 

(see  Figure  l),  the  potential  at  P  due  to  the  electronic  charge  distri- 
bution is  given  by 


27  r' 


I0(r") 


dr 


»t  _ 


N^-   '° 


-  dr' 


y    -12  J  ^12 

This  integral  is  eval\xated  in  Appendix  1,  and  by  proper  substitution  we 
get 


)  ^12 


)1 


dr"  =  e 


(^o  ^  ^) 


1  N  -2r„(ri+a) 


r^M  ■ 


The  potential  at  P,  due  to  both  kinds  of  charges,  is  then 

-  |0(r-)|^  -.   e     e   .  ./.   .  1  ^ /^r^^^  V^^ 


V„  =  —  -  e 
^   ^1 


r^2         r^   V^    V°   ^1-^/ 


The  first  two  terms  on  the  right-hand  side  represent  the  potential  due  to 
two  opposite  point  charges  a  distance  d  apart.  Substituting  |d|  =  a/r^ 
in  this  formula  we  get 


V    =2-  . 

e 

+  e 

A-    1     ^ 

a 

'°       r.- 

-^'of'i  *  h  > 


It  is  seen  that  the  first  two  terms  represent  the  potential  field  due  to 
two  oppositely  charged  point  particles  a  distance  d  apart.  The  third 
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term  is  then  the  effect  of  spead  of  the  electronic  charge  (Figizre  l) . 
Due  to  the  exponential  factor  in  the  third  term,  the  magnitude  of  this 
term,  when  r^  is  either  small  or  large,  is  much  smaller  thaxi  the  difference 
between  the  first  two  terms.  If  we  neglect  this  term  the  above  equation 
reduces  to 

(2.8)       Vp=f- 


'^    V%      ^1(^1-^ 


From  this  equation  we  see  that  for  small  r  ,  Vp  =  Vc  =  e/r^  and  for 

large  r, ,  V  =  2V,  =  ae/r^ .  The  fact  that  V  is  twice  the  value  of  V, 
*   1'  p     d     '  1  p  d 

instead  of  being  equal  to  it  is  due  to  the  way  in  which  we  have  derived 
Vp.  In  this  derivation  we  have  assumed  that  for  large  r, ,  Vp  is  the 
dipole  potential  due  to  a  permanent  dipole  whose  moment  happens  to  be 
equal  to  P  =  e  a/r  ,  as  can  be  seen  from  (2.8).  Actually  in  the  electric 
field  of  an  incident  electron  the  hydrogen  atom  behaves  as  an  induced 
dipole.  A  more  suitable  way  to  derive  (2.8)  would  be  to  consider  the 
Coulomb  and  the  induced  dipole  forces  exerted  on  the  incident  electron, 
and  to  find  the  potential  by  calculating  the  work  done  by  these  forces 
on  the  incident  electron  when  this  electron  moves  from  infinity  to  the 
point  r  .  However, such  a  calculation  is  lengthy  and  is  imduly  refined 
for  the  crude  model  from  which  (2.8)  is  derived.  We  are  content  to 
write  (2.8)  in  the  form 

(2.9)      V„  =  °^ 


1   r^(2r5+  a) 

It  is  seen  that  V  reduces  to  (2.6)  and  (2.7)  for  proper  ranges  of  r  . 
Now, by  assuming  that  the  perturbed  hydrogen  atom  wave  function  is 
given  by  0(r)  =  Ne"'^,  the  potential  field  of  this  atom  can  be  calculated 
and  is  given  by  (Appendix  l) 
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(2.10) 


1   -^''^1 
V0  =  e(,  .  i-)e 


With  two  forms  of  the  potential  of  the  Incident  electron  given  by  (2.9) 
and  (2.10)  we  require  next  that  the  potential  energy  of  the  incident 
electron  in  its  range  be  the  same,  using  either  (2.9)  or  (2.10).  In  this 
way  a  value  for  y   is  obtained.  The  main  reason  for  adjusting  y   and  using 
(2.10)  instead  of  (2.9)  in  the  equation  of  motion  of  the  incident  electron 
is  to  include,  through  7,  the  effect  of  polarization  of  the  hydrogen  atom, 
in  the  exchange  potential,  since  y   appears  in  the  expression  for  this 
potential.  It  is  this  potential,  that  gives  rise  to  a  cross  section  resonance 
when  the  incident  and  the  bound  electrons  have  parallel  spins.  This  is 
discTossed  in  Section  ^. 

If  F  (r..)  represent  the  space  function  of  the  incident  electron,  with 
the  Tipper  and  lower  signs  for  the  total  symmetric  and  antisymmetric  space 
functions  respectively,  the  criterion  that  the  potential  energy  of  the 
incident  electron  for  the  two  potential  fields  be  the  same  is  given  by 


(2.11) 


rj_(2r^-Hi)       ^1 


F"(r^) 


dr^  =  0 


Substituting  the  right-hand  side  of  (2.2)  into  (2.11),  and  noticing  that 


P.  (^i)  are  orthogonel  functions,  we  get 


(2.12) 


je=o 


r^(2r^-KL) 


1  -^'^\ 
1 


i  ('1' 


dTj^  =  0 


This  equation  is  satisfied  if  we  require  that 
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CO 


(2.13) 


r^(2r^-Ki) 


1   "^''^1 
(7   +  7-)e    ^ 


";(^.) 


dr^  =  0  , 


^=  0,1,2, 


The  integration  of  (2.13)  would  give  y   in  terms  of  the  polarizability  a. 

It  is  seen  that  for  each  spin  and  angular  momentum  state  of  the  incident 

+ 
electron,  y   has  a  different  value  and,  furthermore,  since  lL*(r)  is  only 

known  numericeilly,  the  solution  of  (2.15)  requires  numerical  integration. 

A  close  inspection  of  (2.13),  however,  shows  that  y   cannot  be  much  different 

for  different  states  of  the  incident  electron.  The  reason  for  this  is 

that  when  r^  is  very  small  the  sum  of  the  terms  in  the  square  bracket 

approaches  a  constant  and  U  (r, )  behaves  like  r^   (see  Appendix  III.) 

The  contribution  to  the  integral  from  this  region  is  therefore  negligible. 

+  + 

When  U.  (r^)  crosses  the  r^  axis  for  the  first  time  in  a  IL(r^)  versus  r^ 

diagram,  r^  is  of  the  order  of  n/2k  where  k  is  the  wave  nxmiber  of  the 
1  o       o 

incident  electron.  For  small  values  of  k  ,  which  are  of  interest  in  this 

o 

paper,  r^  is  large  enough  to  make  each  term  in  the  square  bracket  in  (2.15) 

sufficiently  small.  We  conclude  that  the  largest  contribution  from  each 

of  the  two  terms  in  the  bracket  to  the  integral  (2.15)  is  from  a  region 

+ 
in  which  ur(r, )  does  not  change  rapidly.  We  therefore  make  the  ajjproximation 

+     ^ 
that  U"(r  )  is  constant  in  this  range  and  write  (2.15)  in  tli®  form 

(2.1^^) 


CO 


a 


p^(2r^  -to) 


'   1   ■^'^l 
(r  +  ^)e   ^ 


dr  =  0 


Integral  (2.1^)  is  evaluated  in  Appendix  II.  The  following  relationship 
for  7  is  found: 

(2.15)       7  =  0.585/a^'^^  . 
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In  Figure  2,  V  and  V^,  as  given  by  (2-9)  and  (2.10),  are  plotted 
q      p 

against  r  for  a  =  2.55  a  and  its  corresponding  vulue  of  7  =  0.44/a 
calculated  from  (2.15).  It  is  seen  that  V  and  Vj  have  the  same  value  at 
the  origin.  As  r  increases  V  becomes  larger  than  V-j  for  the  saone  r.  At 
approximately  r  =  I^Vj  crosses  over  and  becomes  larger  than  V  •  Since  Vj 
vanishes  exponentially,  V  becomes  larger  once  more  for  larger  r  and 
remains  so  for  the  rest  of  the  range. 

2.5  The  choice  of  7  for  different  angulsur  mcanentum  states  of  the  Incident 
electron 


By  means  of  (2.15)  the  veilue  of  7  depends  on  the  polarizability  of  the 
hydrogen  atom.  So  fax  no  precise  measurement  of  this  polarizability  has 
been  made.  Even  a  precise  measurement  of  the  polarizability  would  not 
give  a  satisfactory  vaJ-ue  for  7.  The  reason  for  this  is  that  at  the 
classical  distance  of  closest  approach  of  the  incident  electron  the  electric 
field  in  the  vicinity  of  the  hydrogen  atom  is  highly  inhomogeneous  and  its 
magnitude  is  far  beyond  the  magnitude  of  the  electric  field  used  to 
measure  polarizability  of  this  atom.  The  maximimi  uniform  electric  field 
used  to  measure  this  polarizability  has  been  500  e.s.u.  ^  -'  vhile  the 
incident  electron  at  a  distance  equal  to  a  Bohr  radius  from  the  nuclevis 

establishes  an  electric  field  at  the  nucleus  equal  to 

-10  -Q  ?  7 

h.Q   X  10   7(5.29  X  10  ^)   =  1.7  X  10'  e.s.u.  Since  the  formula  p  =  oE 

is  a  linear  approximation,  this  method  of  determining  7  is  only  an 

approximate  one.  It  is  evident,  however,  that  7  depends  on  the  amount  of 

penetration  of  the  incident  electron  in  the  potential  field  of  the  atom 

and  is  therefore  related  to  the  classical  distance  of  closest  approach. 


Electrostatic   potential  in  atomic  units 
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This  di stance  is  given  by 

(2.16)        d  =  L/p  =  /JWTJ  ix/ik^  =  i/k^, 

L  and  p  being  the  angular  and  3j.near  nKmentum  of  the  electron.  While  for 
small  >f  this  distance  is  of  the  order  of  the  Bohr  radltis,  its  value  is 
different  for  different  ><^.  Hence  7  should  be  different  for  different  angular 
momentum  quantum  numbers.  The  existence  of  the  electrostatic  exchange 
potentieuL  as  derived  in  the  next  section  makes  the  value  of  7  different  for 
incident  electrons  with  the  same  angular  momentinn  but  different  spin  states. 
For  a  singlet  state  vhere  the  two  electrons  have  antiparallel  spins,  the 
exchange  potential  turns  out  to  be  a  repulsive  one.  For  a  triplet  state 
where  the  two  electrons  have  parallel  spins  the  opposite  situation  takes 
place:  the  two  electrons  have  a  small  probability  of  being  close  to  one 
another,  and  this  gives  rise  to  a  total  potential  more  attractive  than  the 
direct  potential  alone.  In  this  way  the  exchange  potential  modifies  the 
total  effective  potential  that  the  incoming  electron  sees,  and  this  in 
turn  changes  the  distance  of  closest  approach  of  the  two  electrons  with 
a  corresponding  change  in  7.  In  the  same  way  it  can  be  inferred  from 
(2.16)  that  7  depends  on  k  .  We  conclude  that  7  covers  a  range  of 
veJ-ues  for  different  states  of  the  Incident  electron.  In  this  paper  7  is 
fixed  as  follows: 

As  was  discussed  earlier  in  this  section,  the  partial  wave  for  the 
s-scattering  is  spherically  symmetric  about  the  origin,  so  that  on  the 
average  the  atom  is  not  polarized.  We  can  thus  take  7=7  =  l/a  for 
8-8cattering.  With  7  fixed  for  the  s-state,  it  is  chosen  for  the  other 
angular  momentum  states  so  that  it  leads  to  a  value  of  the  total  cross 
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section  in  closest  agreement  with  the  experimental  results.   In  this  way 
we  find  that  y  =   0A63/a  .  We  anticipate  the  results  of  the  numerical 
integration  and  mention  that  with  this  choice  of  y   the  cross  section,  due 
to  the  states  higher  than  Jl  =  1  ,   is  negligible.  By  (2.I5)  7  =  0.^63/a 
corresponds  to  a  =  2.00  or    =  2.96  X  10  ^cm  according  to  our  model. 

Sheffers  and  Stark"-  -*  have  mesisured  the  polarizability  of  the  hydrogen  atom 

-25  5 
and  have  found  a  value  equal  to  (5  +  l)  X  10   cm  for  an  electric  field 

intensity  of  7  X  10  V/cm.  It  is  seen  that  the  value  used  in  this  paper 

agrees  with  the  experimentaJ.  value  within  the  limits  of  the  experimental 

accuracy.  The  theoretical  value  for  the  polarizability,  however,  is  higher. 

3  -25   3  fsl 

Its  value,  based  on  perturbation  theory,  is  4. 50  a  or  6.77  X  10    cm  '-  -"  . 

Further  discussion  of  the  value  of  y   is  given  in  Section  5- 

5.  The  wave  equation  of  the  incident  electron 

With  a  known  form  of  one  of  the  product  wave  functions,  it  is  possible 

to  obtEdn  the  wave  equation  of  the  incident  electron.  The  total 

Hamiltonian  of  the  system  of  two  electrons  and  one  proton  is 

^2  ^   ^2  _    2    2    2 
2m  1   2m  2   r,   r^   r.p 

where  m  and  e  are  electronic  mass  and  charge,  respectively,  r^  and  r^  are 
position  vectors  of  the  incident  electron  and  bound  electron,  respectively, 
r^   =  jr  -  r  I .  E  is  the  total  energy  of  the  system.  The  total  wave  function 
of  the  two  electrons  including  their  spin  coordinates  should  be  anti- 
symmetric with  respect  to  the  interchange  of  these  electrons  in  order 
that  this  wave  function  satisfy  the  Pauli  exclvision  principle.  This 
condition  leads  to  the  symmetric  space  function  for  the  singlet  spin  state 
in  which  the  electrons  have  antiparallel  spins  and  to  the  antisymmetric 
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space  function  for  the  triplet  spin  state  in  which  the  electrons  have 
parallel  spins.  The  space  fimctions  are  then  of  the  form 


t,^    ^ 


(5.2)    4''(rj_,r2)  =  F"(r^)0(r2)+  F'(r2)0(r^) 

•where  the  upper  sign  corresponds  to  the  symmetric  and  the  lower  sign  to 

the  antisymmetric  space  function,  respectively.  We  introduce  a  new 

I  2  2 

constant  7  =  1/a  =  me  /fi  where  a  is  the  cleissical  Bohr  radius  and  let 


(5.1)  operate  on  (5«2),  obtaining 


(5-3) 


^+V^  +  27(—  +^-  i— )  + 
1    2    'o^r^   r^  T^' 


2m 


1    2    ^o^r^   r^   r^^   ^2 


E 


F~(r^)0(r2) 
F"(r2)0(r^)  =  0. 


The  wave  equation  of  the  incident  electron  is  obtained  by  multiplying  (5«5) 


tiy  0(rp)  and  integrating  over  all  coordinates  of  r^, 


i3M 


f 


J 


0(r2) 


1    2    ^o^r^   r^   r^^ 


2m 


F'(r^)0(r2)d?. 


1   0(^2) 


1    2    'o^r,   r^   r,„*'   ^2 


E 


12    * 


F"(r2)0(r^)dr2  =  0 


The  integration  is  carried  out  by  noting  that  0(r)  =  Ne"'^.  The 
normalization  constant  N  =  27^'  /  \J^     is  obtained  by  requiring  that 
I  o)2(j.  )(j3.  =  1.  Then  the  following  relationships  hold: 


r 


(5.5) 


[- 


27, 


Ci(r)  = 


"2  ^^yy)" 
7  +  — = — 


0(r) 


(5.6) 


0(r) 


r 


0(r)dr  =  27^7  -  j" 
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Making  use  of  (5'5)  and  (5.6)  in  (5.^)  we  now  have 


(5.7) 


^.2r,r-r^.pE.2r  f|0(r,)|^(i-.i-)d^^ 


0(r2) 


^^ 


F-(r^) 


F'{r^)0{r^)6r^ 


Letting  k  ,p  represent  the  wave  number  and  the  momentum  of  the 
incident  electron  at  large  distances  from  the  atom,  respectively,  and 

noting  that  the  energy  of  the  normal  state  of  the  hydrogen  atom  is 

2         2 
-e  /2a  =  -e  7  /2,  we  obtain 
'   o      'o'  ' 

(5.8)    k^  =  p^/fi^  =  2m(E  +  e^yj2)/ii^   =  2mE/h^  +  Yq  ' 

where  E  is  the  total  energy  of  the  system.  Also,  by  Green's  second  theorem, 

(3-9)     f  \0(.T^)'^-(i^)  -   p-(l?2)^0(r2)]d?2 

The  second  integral  in  (5.9)  is  carried  over  the  stirface  of  a  large  sphere 
of  radius  r-  ->  00,  and  since  0(r)  is  botmded,  this  function  and  its 
derivative  vanish  on  the  sTorface  of  this  sphere.  The  surface  integral  in 
(5 '9)  then  vanishes  and  we  get 

(5.10)     J  0(r2)^"(r2)<3^2  =  j  ^'(^^)'^0ir^)<^2* 

Making  use  of  (5.8),  (3. 10),  (5.5)  and  the  resiilt  of  Appendix  I  we  finally 

obtain  from  (5-7): 


(3.11) 


■27r, 


.1^  ^o-  ^ V  ''^'■^  2rj7  -  ^)e         \ 


l^Cr^)  0(r2) 


"2    P  2.  ^(yy^ 

^0^27  -7^+  —^ 


27^  +,_v  ,  ^ 


r{t^)^^mr^)j   0(t^)j^  r{r^)dT^ 


=   0 
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Eqoiation  (5.II),  except  for  the  middle  term  and  the  third  term  in  the 
bracket  of  the  first  term,  is  the  Eartree  Fock  equation  for  the  system  of 
two  electrons  moving  in  the  field  of  a  nucleiis  •-  -• .  These  terms  vanish 
if  0(r)  represents  the  groimd  state  of  the  hydrogen  atom  and  if  the 
wave  fimctions  of  the  two  electrons  can  be  made  orthogonal.  The  former 
equation  was  derived  without  minimizing  the  total  energy  of  the  system 
with  respect  to  the  two  wave  functions  in  pontrast  to  the  usual  derivation 
of  the  Hartree  Fock  eq-uations. 

The  three-dimensional  equation  (j.ll)  can  be  reduced  to  a  one -dimensional 

eqioation  by  expanding  F  (r)  in  terms  of  the  zonal  hannonics, 

+  _,    ,  00  + 
(2.2)   F"(r)  =^Y1  u:(r)P  (n)  ; 

l/r  p  can  be  expanded  in  a  similar  way: 
^         ^  CO  Tg  n 
12    1  n=o  1 


^J>'l2l 


m=l  ^ 


where  [x       =  cos  0  ,  n^  =  cos  0^,  p.^   =  cos  9^,   0^^  ^^i^^g  *^e  angle 
between  r^  and  r^,  9^,  0^  and  6^,   0^  are  the  spherical  coordinates  of  r 
and  r  ,   respectively.  With  the  help  of  these  expressions  the  following 
equations  result: 
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-2rr, 


V?  +  27  (7  +  — )e         ^  +  k^  -   (7  _v) 
1  o  r  '  00' 


F"(r^) 


00 

H 
£=0 


-27r. 


^  +  27j7+^)e         h^l-   (7^-7) 


7^  u;(r^)P^  (,^) 


1     ^ 
1  Ko         -^ 


";(^i) ' 


^12  2 


OD        00 

j2  =0  n=o 


/  ^1    /  ^1    /  2n 


o      '-1     'o  ^  '1 


y 


.00        +1    ,  2rt 


0(rj 


u;(r^) 


2'        r, 


-1    ''o 


2  r: 


^  P/m2)V^i)V^2) 


r2dr2(iM2'i02 


00 


/'l 


0(r   )u"(r   )r  -^'^■'■dr     +  r 
"^     2    i>  ^    2     2  2         1 


00 
^'1 


OD 


0iT^)F'ir^)dT^  =  i^n  0(r2)U^(r2)r2(ir2     , 


f  0(r  )F"(r  )dr  / 
§ 2 2   ^  i^^ 

"2  J, 


00 


0{T^)V'^{T^)dT^    . 


If  we  substitute  the  right-hand  sides  of  the  last  four  equations  into 
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(5-11)  and  equate  the  coefficients  of  P-yCu-,)  to  zero,  the  following  equations 
result  : 


r. 


+  2rj7  +  -)e    -  -^  +  k^  -  ( vr) 

dr^  1  r^ 


u;(r,) 


+  4n5   0(r, ) 


00  ,  ^  J. 

(k^+  27^-  7^)  J   <l>{v^)\i'{T^)T^^^M7^-7)\        <l>{v^)\i^{r^)^^ 


^7, 


2i+l 


I  «'i' 


r-^i 


-e*i 


«((r,)u;(r3)rf  V  "•/ 


/  oo       _^      ^^ 


=  0 


.5/2 


^=  0,1,2,..., 


Since  p(r)  =  27-^'  /  ^57  e"''^,  these  equations  can  be  vn^itten: 


(5.12) 


dr. 


^     f       1   ^    "^''^1      J^(^+l)   ^  ^2       /  n2 

+  27  (7+  — )e  ^    „       +  k     -   (7  -7) 

or,  2  o         'o 


u;(r,) 


+  J+7^6     r,e       ^ 


00  ,00 

(k%27''-7^)  e       '^^(r2)r2dr2+^(7^-7)        e       'lf{T^)dx^^ 


87^     -7r^ 
+  2J?+1     ^ 


r,      j?+l 


/^    1  r^  ,^2 


7T^  + 


e       ^U^(r2)dr2+ 


X  CD        jg+l 

f         r,  -7r„  + 


V"        '^i("2)^2 


=  0 


>t     ~      Wy±,,C,.c.. 


The  factor  (7  -7)  appearing  in  the  first  term  of  (3.12)  corresponds  to 
the  amount  of  energy  that  the  incident  electron  has  spent  in  polarizing 
the  hydrogen  atom.   It  owes  its  existence  to  the  fact  that  we  have  not 
properly  included  the  boundary  conditions  in  our  problem,  namely  that  when 


r  approaches  infinity  the  wave  function  must  reduce  to  an  outgoing 
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spherical  wave  times  the  ground  state  hydrogenic  wave  function.  Presumably, 

if  the  boundary  condition  were  properly  taken  into  account,  such  a  change 

of  the  energy  scale  would  not  occur. 

Nevertheless,  some  decision  has  to  be  reached  about  what  to  do  about 

this  term^ since  for  low  energies  of  the  incident  electron  and  for  7*7, 

we  do  not  obtain  a  solution  which  is  outgoing  at  infinity.  For  s-scattering 

we  assume  that  7=7  and  the  problem  disappears.  For  higher  angular 

momentum  states,  when  7  ^  7  ,  we  can  see  that  effectively  the  energy  of  the 

incoming  particle  has  been  reduced  to  a  value  equal  to  |k  -  (7-7  )  1  ' 

instead  of  k  .  This  is  correct,  the  difference  being  accoimted  for  by  the 

fact  that  some  of  the  energy  of  the  incoming  particle  has  been  used  to 

polarize  the  atom.   If  we  had  a  proper  model  for  depolarizing  the  atom  as 

the  incoming  particle  recedes  we  would  then  have  the  energy  reduced  again 

to  k  .  Our  feeling  then  is  that  the  term  in  7-7  is  not  one  which  would 
o  o 

persist  in  a  correct  model  and  thus  we  have  eliminated  it  from  the 
equation.  At  worst  the  effect  of  what  we  have  done  is  to  somewhat  distort 
the  energy  scale  of  the  incoming  particle. 

h.     Numerical  integration  of  the  radial  equation 

^.1  Method  of  solution 

To  facilitate  the  numerical  integration  we  introduce  the  following 
dimensionless  variables: 

(U.l)    X  =  k^r^,   x-  .  k^r^,   P„  -•  yj\  ,      P  .  yl\   . 

Equations(5.12),  with  the  omission  of  the  factor  (7^-7)  in  the  first  term 
and  some  rearrangements  of  the  limits  of  inegration,  can  now  be  written: 


22  - 


(^.2) 


dx  X 


v;m 


2i+l 


r  X 


,x+i 


£ 


.^+1 


■^^'u;(x  )dx  -  c^-x^^i 


=  0 


(^•5) 


00 


^i-= 


,-je 


5 
2P„ 


(1+2P^-P^)x'+  i^(p^-P)  . 


e"^^  l^"(x')dx' 


/  =  0,1,2, ... 
+ 
Eqiiatlon  (^.5)  shows  that  the  values  of  C.  cannot  be  determined  lonless  the 
+  *  + 

ur(x)  are  known,  and  (^.2)  shows  that  the  differential  equation  for  Urt(x) 

cannot  be  integrated  unless  the  C^  are  known.  One  way  of  overcoming  th:.s 

difficulty  is  to  assume  some  suitable  values  for  C. ;   with  these  values  (^.2) 

+ 
and  (^.5)  are  integrated,  and  with  the  new  values  of  the  C.  the  same 

+ 
operation  is  repeated  until  consistent  values  for  the  C.   are  obtained.  A 


more  siiitable  method  is  as  follows: 

+ 


We  express  U"(xJ  in  terms  of  the  two  other  functions  U  (x)  and  IL  (x.l 


+ 

i. 


such  that 


+  + 


(k.k)         V'Ui    =  Vj   U)   +  C"  V'   (x)   , 

+         + 
where  U'  ^x)  and  U7  (x)  satisfy  the  following  equations; 


(^.5) 


dx" 


(p,  l)e-P^-  ^^  .   1 


\^^' 


(^ 


x 


i!+l 


r^""  u"  (x')dx' 

-*1 


=  0 


and 
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i^.6) 


7  2  ^  2P, 

cLx 


(p  .  l)e-P-.  A(.^  ,  1 


U"  (x) 


2J?+1 


>o 


,/+l 


^+1 


x^     '£ 

X 


If  we  miatiply  ('+.6)  by  C"  and  add  to  (^.5),  we  obtain  (^.2). 

The  value  of  C^  is  obtained  by  substituting  C^-^)  into  (^.5)  which 
gives 


oo 


i^-l) 


x'   -  ^|(l+2p2-P^)x'+  i^O^.p) 


e"^  U"  (x  )dx 
■'  1 


r 


00 


,-i  8 


2P, 


(l+23^-3^)x'+  '+(P^-P) 


e"^""  u:  (x')dx' 
■^2 


+      +  + 

For  large  x,  U  (x),  U.  (x)  and  U  (x)  are  sinusoidal,  each  with  its 

own  amplitude  and  phase  shift.  The  asymptotic  form  of  (^.^)  is  then 


+  + 


(h.Q)  A'sinfxf  6")  =  a"  sin(xf  6~  )  +  C"a"  sin(xr  6"  ) 

When  this  equation  is  expanded  we  have  the  identities 


+  + 


A"sin  5"  =  a:  sin  6"  +  C'aT  sin  5 


A^cos  5^  =  A^  cos  5^  +  C^aJ^cos  5^ 


+    +  + 
1   ^  ^2 


or 


(i^.9)   tan  5^  = 


+  +  +  +  + 
a"  sin  5.  +  C'aJ  sin  5_, 

+  +  +  +  + 
a"  cos  67  +  qTa"  cos  8" 
^1    /l   ^i'2    ^2 


Equations  (^-5)  and  {h.G)   can  be  integrated  from  the  origin  with 


2k 


+ 


proper  boimdary  values.  C~  is  then  calcxolated  by  means  of  (^.7).  Using 

■^   +         + 
the  asymptotic  anrp3j.tudes  of  U~  (x)  and  U^  (x)  together  with  their  phase 

+   -''1        ^2  + 

shifts  and  the  values  of  c7,  we  obtain  the  phase  shifts  of  ur(x)  from 

('^-.9)-  To  find  the  total  elastic  scattering  cross  section  we  note  that 

5/^  of  all  the  incident  electrons  have  spins  parallel  to  the  bound 

electron  spin  and  the  remaining  l/k   have  antiparallel  spins  ^     ^  .  The 

ri2l 
total  elastic  cross  section  is  then  given  by*-  -• 


.   oo 

(4.10)   S  =^XI  (2i+l) 

k"^  i-o 
o 


1.2^+   5.2- 
j^   sin  5^  +  ^  sm  ^ 


In  this  paper  the  numerical  integration  ha^  been  carried  out  for  j?=  0,1,2. 

+ 
The  integration  cannot  be  started  without  the  values  'of  UJ  (x)  and 
+  ^1 

U/,  (x)  at  the  origin.  These  have  been  calculated  in  Appendix  III.  The 

results  are 

(4.11)      u;  (x)  =    b;x^-*-i 

^  1   x->o  ^  1 


+ 


iop  p5,<+5 


+  + 

The  constant  bZ   in  (4.11)  is  arbitrary.  This  makes  the  ajnplitude  of  UT 

and,  through  this,  the  amplitude  of  LL  arbitrary  at  large  distances.  This 

is  unimportant  since  we  are  concerned  with  the  values  of  the  phase  shifts 

+  + 

of  U-  only.   It  should  be  noted  that  the  amplitude  of  U~  is  fixed  by  (4.12) 

*  ActuaJJLy  the  integrals  appearing  in  (4.7)  are  evaluated  during  the 
coiirse  of  integration  of  (4.5)  and  (4.6). 


h . 2     Method  of  intogration 

•  ri3i 

We  use   Millies L  -^  method  to  integrate  the  equations.  If  a  second-order 
differential  equation  of  the  following  form  is  given 

(i+.l6)    y"  =  f(x,y) 

the  first  fouT  points  are  integrated  by  means  of  the  Taylor's  series  to 
initiate  the  solution.  Then  the  formula 


(^•1^)     ^n.l-  2^1^  ^n-5  =  '^^'^Cl^  5  ^'y;.,)-^  id  ^^y^°^ 

2 
is  used  to  find  y  , -5  y"  in  this  formula  is  defined  as 
n+1 

C^.is)      s'y^  -  yn.i  -  2yn -^  Ci  • 

Since  the  values  of  y  at  the  four  points  x  ,,  x  „,  x  ^  and  x  are  knovm, 

n-^  n~d       n-1      n 

(4.17)  can  be  used  to  find  y  , .  As  a  check  on  the  accuracy  of  this  value 
the  following  correction  formula  is  used: 

2      1   2       h^(^) 
(^.19)    y„,,-  2y^  -  y„.,  =  h2(y^  .  i^  5\"  )  -  -^     , 

to  find  a  more  correct  value  of  y   .   If  the  difference  between  the  two 
values  is  less  than  a  prescribed  value  for  the  error  of  integration,  the 
integration  can  proceed  to  the  next  point. (The  interval  of  integration 
can  even  be  doubled  to  accelerate  the  integration.)  But  if  the  difference 
is  more  than  a  prescribed  value  the  interval  is  halved  and  the  integration 
is  repeated  until  the  desired  accuracy  is  attained.  From  (4.17)  and  (4.19) 
this  difference  is  17h°  •  y^  V24o. 

The  integrals  appearing  in  (4.5)  and  (4.6)  can  be  integratec  from  the 
origin  while  the  integration  of  the  equations  is  progressing.  Simpson's 
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X 

/  2 


y'(x)(ix,  then 

X 


rule"-  -•  is  used  for  evaltiat  ion  of  the  integrals.  If  y  = 
according  to  this  rule 

h  h5v(5) 

(^.20)   y^-  y^  =  f(y;  +  ky'^  ^  y')  -  ^^     . 

This  rule  is  not  only  accurate  but  simple  because  a  knowledge  of  only  one 
previous  value  of  the  integral  is  reqxiired  to  continue  the  integration. 

Equations  (^.5)  and  (^.6)  were  integrated  on  the  IBM  704  digital 
conrputerv 

5.  Interpretation  of  the  results  and  ccmparison  with  experiment 

As  v&B   pointed  out  in  Section  2,  the  choice  of  7  =  l.OO/a  for 
s-scattering  and  7  =  0.465/a  for  higher  angular  momentum  scattering  wovild 
lead  to  a  total  elastic  cross  section  in  close  agreement  with  the 
experiment.  The  reason  for  taMng  7  =  l.OO/a  for  s-scattering  is  that  the 
s-partial  wave  is  spherically  synnnetric  about  the  origin,  and  polarization 
of  the  atom  can  be  neglected.  There  are  other  reasons  for  this  assumption. 
The  s-scattering  cross  section  at  k  =  0  increases  with  a  decrease  in  7 
according  to  the  formula 

(5.1)     S(k^=  0)  =  l^na^^l/7^  , 

[15! 
S  being  the  s-scattering  cross  section  and  a  the  scattering  length"-  -* . 

For  intermediate  values  of  k  the  largest  contribution  to  the  phase  shift 

is  from  the  Coulomb  potential^ and  polarization  does  not  play  an  important 

role  while  in  other  angular  momentum  states  the  Coulomb  potential  is 

effectively  screened  by  the  "centrifugal  potential*  and  the  effect  of  the 

[15]* 
polarization  potential  is  more  evident'-  -■  .  These  considerations  show 


*  Similar  conclusions  are  reached  with  the  model  used  in  [}] ,   Section  5- 
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that  the  only  importaxit  effect  of  taJcing  y   less  than  1  in  the  s-state  is 

to  increase  the  total  cross  section  at  k  =0.  Since  Bederson's 

o 

experimental  curve  does  not  extend  to  k  =0,  we  cannot  determine  y   in  this 

case  by  comparing  our  results  with  the  experimental  ones.  However,  the 

calculated  cross  section  according  to  the  present  model  with  y   =  l/a 

agrees  with  the  independent  calculations  of  the  elastic  scattering  cross 

section  of  electrons  by  hydrogen  atoms  at  k  =0^  ■^''-  -' .  We  have  therefore 

adopted  the  value  of  7  =  7  for  s-scattering  in  this  paper.   For  other 

angular  momentum  states  it  is  found  that  7  =  0.^63 /a  gives  rise  to  a 

curve  for  the  total  cross  section  in  close  agreement  with  the  experimental 

curve  of  B.  Bederson. 

With  this  value  for  7  the  angular  momentum  states  higher  than  i'  =  1  do 

not  contribute  appreciably  to  the  cross  section.  For  exaaiple,  the  d-state 

has  a  maximum  cross  section  about  2jta  when  k  =  O.T/a  •  For  this  reason 

o      o     '  o 

the  contribution  of  angvilar  momentvmi  states  higher  than  /=  1  are  neglected. 

In  the  following  pages  the  resxilts  of  the  nimierical  integration  are 
tabulated  and  are  shown,  graphically. 

In  Table  1  are  the  numerical  values  of  the  phase  shifts  and  the 
scattering  cross  sections  when  the  incident  electron  is  in  each  possible 
spin  and  angular  momentum  state,  and  for  7  =  1.00.  In  Table  II  the  same 
results  are  given  for  7  =  O.it-65  for  p- states. 

Figure  \   shows  the  phase  shifts  as  a  function  of  the  incident  electron 
energy  for  the  2  p-states  considered  in  this  paper.  It  is  seen  that  the 
singlet  phase  shift  assumes  negative  values  for  some  range  of  the  incident 
electron  energy.  This  is  explained  by  the  fact  that  in  the  singlet  state 
the  two  electrons  have  some  affinity  for  each  other.  This  affinity  is 
expressed  by  a  positive  electrostatic  exchange  potential  energy  as  is  seen 
by  the  lower  sign  in  (5. 11).  This  potential  energy  may  be  large  enough 
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ELASTIC    SCATTERING    OF   ELECTRONS 
BY     ATOMIC    HYDROGEN 


Phase  shift     of  the 
s-scattering  when  y  -  I 

I.  Phase  shift  of  the  triplet  state 
as   obtained  by  Morse  8    Allis  and 
in  this  paper. 

Phase  shift  of  the  singlet  state  as 
obtained  by  Morse  8  Allis. 
Phase  shift  of  the  singlet  state 
as  obtained  in  this  paper. 
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FIGURE    5.    ELASTIC  SCATTERING   OF  ELECTRONS 

BY  ATOMIC    HYDROGEN 


S.  n!   Cross  section  of  the  incident 


electron   in  X  spin  state    and 
£angular    momentum  state. 
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to  outweigh  the  attractive  direct  potential  energy.  The  zero  phase  shift 
signifies  that  the  average  exchange  potential  completely  balances  the 
average  direct  potential. 

The  phase  shifts  for  the  triplet  states,  on  the  other  hand,  are  always 
positive  due  to  the  fact  that  the  exchange  potential  in  these  states  is 
always  attractive  and  adds  to  the  effect  of  the  direct  potential. 

In  Figure  5  the  contribution  to  the  total  cross  section  from  each  state 
is  shown  graphically.   It  is  seen  that  the  largest  contribution  comes  from 
the  triplet  state  with  Jl  =  \. 

In  Figure  6  the  total  cross  section  is  shown  for  7  =  l.CX)  in  the  s  and 
p-states  (cvirve  l);  and  for  7  =  1.00  in  the  s-state  and  7  =  0.^65  in  the 
p-state  (curve  2).  The  experimental  curve  of  B.  Bederson  is  8Q.S0  redrawn 
for  comparison.  It  is  seen  that  when  7  changes  from  1.00  to  0.^65  in  the 
p-state,  the  theoretical  curve  reproduces  most  of  the  features  of  the 
experimental  curve.  As  shown  in  Section  2,  a  vedue  of  7  =  0.i|-65  corresponds 
approximately  to  a  value  for  the  polarizability  of  the  hydrogen  atom 
equal  to  2.00  eP   or  2.96  X  lO"  ^  cm^.  Sheffers  and  Stark"-  -■  have  measured 
polarizability  of  the  hydrogen  atom  and  have  found  a  value  equal  to 
(3  -  1)  X  10"^^ cm^  for  an  electric  field  of  7  X  10  v/cm.  It  is  seen  that 
the  value  used  in  this  paper  agrees  with  the  experimental  value  within 
the  limits  of  the  experimental  accuracy. 

Figure  5  is  the  ccorparison  of  the  s-scattering  phase  shift  for  7  =  1.00 
as  obtained  in  this  paper  and  as  calculated  by  Morse  and  Allis  •■  -■  .  There 
are  some  discrepancies  for  singlet  scattering  due  to  the  errors  pointed 
out  in  page  5.  This  discrepsuacy  may  explain  most  of  the  differences 
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for  the  singlet  state  between  the  curve  obtained  by  a  variational  method  by 

ri9i 

Massey  and  Moisewitsch'-  ->   and  the  curve  obtained  by  numerical  integration 
by  Morse  and  Allis. 


It  should  be  noted  that  a  better  agreement  with  the  experimental  curve 
can  be  obtained  if  contribution  to  the  totaJL  cross  section  by  higher 
angixlar  momentum  states  with  different  7  is  taken  into  account.  These 
values  for  7  shoiald  be  less  than  OA63.  As  J   increases,  the  classical 
distance  of  closest  approach  increases  and  the  incident  electrons  with 
lower  values  of  k  do  not  enter  the  potentieil  field  of  the  atom.  As  a 
result  the  maximum  contribution  to  the  cross  section  appears  at  higher 
values  of  k  .  In  this  way  a  theoretical  curve  in  better  agreement  with 
the  experimental  curve  can  be  obtained.  We  have  not  carried  out  this 
program  here  principally  because  there  is  no  unique  way  of  finding  the 
values  of  7  for  each  state. 

We  also  enumerate  the  other  approximations  ixsed  in  the  present  static 
model  which  may  be  partly  responsible  for  the  discrepancy  between  the 
theory  and  the  experiment. 

1.  The  charge  distribution  of  the  bound  electron  depends  on  the  position 
vector  r^  of  the  incident  electron.  In  the  assxaned  model  this  charge 
distribution  is  independent  of  r^  and  is  a  fxmction  of |rp|only. 

2.  While  the  perturbed  hydrogen  atom  wave  function  is  spherically 
symmetric  with  respect  to  the  origin  in  s-scattering,  it  has  a  different 
symmetry  in  higher  angular  mcanentum  scattering.  We  have  assumed  that  this 
wave  function  is  spherically  symmetric  in  all  states. 
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5-  The  neglect  of  the  term  (7  -7)^  in  (3.12)  has  the  effect  of 
distorting  the  energy  scale  in  some  iinknovn  way. 

k.     A  less  important  approximation  was  used  in  the  derivation  of  (2.l4) 
from  (2-11)  by  assuming  that  the  radial  wave  functions  do  not  change  rapidly 
in  the  course  of  integration. 

5.  The  total  wave  function  is  given  only  approximately  as  a  symmetrized 
product  of  one-electron  wave  functions. 

6.   Conclusion 

By  maJcing  some  reasonable  assumptions  in  the  scattering  of  electrons 
by  hydrogen  atoms  it  was  shown  that  the  essential  feat\ires  of  the 
experimental,  resiilts  of  B.  Bederson  can  be  reproduced  as  is  shown  in 
Figure  6,   Section  5-  The  discrepancy  between  the  two  curves  is  thought  to 
be  due  to  the  approximations  made  in  the  assumed  model.  These  approximations 
are  enimierated  at  the  end  of  Section  5'  The  resonance  in  the  elastic  cross 
section  of  electrons  by  hydrogen  atcans  at  about  5  ev  incident  electron 
energy,  as  described  by  B.  Bederson,  is  believed  to  be  due  to  a  triplet 
state.  The  largest  contribution  to  the  resonance  in  this  state  comes  tram 
M.=  \   although  higher  angxilar  momentimi  may  also  be  effective. 

The  model  outlined  in  this  paper  has  not  the  superiority  of  a  model  in 
which  the  bound  electron  wave  function  changes  continuously  with  the 
position  vector  r  of  the  incident  electron  but  this  disadvantage  is 
outweighed  by  its  simplicity.  This  model  caji  be  extended  to  the  scattering 
of  electrons  by  any  atom,  provided  the  potential  field  of  the  atom  is 
defined. 
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Appendix  I 

The  potential  due  to  the  wa.ve  function  0(r)  is  readily  obtained  by 

making  use  of  the  Gauss '  theorem.  The  charge  density  of  the  bound 

2      p  -27X2 
electron  is  -  e|0(r2)|   =  -ell  e     .  The  contribution  to  the  potential  of 

the  incident  electron  from  all  negative  point  charges  in  space  for  which 

r^  >  rp  is  the  same  as  the  contribution  from  all  these  charges  placed  at 

the  center  of  the  coordinate  system.  The  contribution  to  the  potential 

from  all  point  charges  for  which  r^  <  r^  would  not  change  if  the  incident 

electron  were  placed  at  the  origin.  Prom  these  considerations  the 

following  resiilts  are   obtained: 


'^ d:^^  =  _lirteN^(^) 


^12      2         ^T^ 


J 


e      r^dr^-WelT   e     ^"^^2*^2 
o  ''r 


where  in  the  last  equation  N  =  27^'  /  \f^  has  been  substituted;  this  is 

the  result  of  normalization  of  0(r)  to  unity.  Taking  the  potential  due 

to  the  proton  into  account,  the  total  electrostatic  potential  acting  on 

the  incident  electron  due  to  the  pertvirbed  hydrogen  atom  with  the  wave  function 

0(r)  is  given  by 


^p-  =  /%-^)l^(-2)l'4 


1   "^^1 
=  e(7+  — )e 

1 
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Appendix  II 


Evaluation  of  the  integral: 
00  r- 


a 


r(2r''-Hx) 


(7  -  i)e-2''- 


o  "- 
Integrating  the  first  term: 


dr  =  0  . 


00 


r(2r-'-Kx) 


dr  =  a 


1  ,    r^ 


2r-' i-a 


00 


[20] 


=  -  ^  log  2 


(log  r)  +   -  log  a 
r->o  -^ 


Integrating  the  second  term: 

f  " 

-  7 


A 


e  '  dr  =  77^  =  -  - 
2y  2 


Integrating  the  third  term: 
/    -27r 


^       dr  -  C  +  log  27  +  (log  r)^_^ 


[21] 


where  C  =  0-577  is  the  Eioler's   constant.     Adding  the  three  terms 
^  log  I  -  I  +  C  +  log  27   -  0   , 


or 
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Appendix  III 

+      + 


Boxmdary  Values  of  IL      and  IL 

1     ^2 


By  Taylor's  expansion  of  II,  about  the  origin j(^. 5)  takes  the 

1 
following  from  for  small  values  of  x : 


(1) 


^2  23^ 

7^  ^  IT 

dx 


U^   (x)    =  0     , 


/=  0 


(2) 


i(i+l) 


^2  -2 

dx  X 


u:   (x)    =  0     ,         i  f  0   . 

-^1 


OD 


Substituting  UT    (x)    =  ^  a  x     in  (l)  we  get: 
•^1  v=o  ^ 


03 

a     =  0  ,        5~"  a 
o  ' — ,    V 


v=i     L 


v(v-l)x"-2  +  2g^x^-^ 


=  0 


This  indicia!  equation  is  satisfied  if 


a  ., •  v(v-l)  +  2P  a  =  0  , 
v+1   ^         o  V 


whence 


=  '^°   a 
v+1    v( V-1)   V 


Then  all  the  a  are  finite,  and 

(5)      U"   (x)  =  a^x  . 
"l   X->o  - 


Substitution  of  the  same  expansion  into  (2)  resiilts  in 


00 


a     =  a, 
o 


^1  =  0,     YL 

with  V  =  J?+  1  and 


a 

v=2   '^ 


v-2 


vCu-Dx^-'^-iCZ+Dx 


v-2 


=  0 


'"     "/V='=c-V^-^ 


jM 


3^ 


It  is  seen  that  (h)   includes  (3)  as  well  and  then  is  valid  for  all  values 

of  JP  . 

+ 

To  find  the  boimdary  value  of  U^  we  use  again  the  Taylor  expansion 

+  2 

of  U^  in  {k.6).     For  small  values  of  x  this  equation  reduces  to 


(5) 


,2    2P   + 
d      o 

+ 


ax 


2    X 


U"  (x)  - 


21+1 


.M 


=  0  ,    i=  0 


(6) 


d^     i(i+l) 

2  '    2 


+  83  B^    g 


The  homogeneous  solution  of  (5)    end  (6)   are  the  same  as   (l)   and  (2).     V/e 

only  need  the  inhomogeneous  solution  of  (5)   and  (6)   to  satisfy  the 

+  +00 

boundary  values   of  tU.     Substituting  the  expansion  U^   (x)    =  ^  a  x 

-^2  v-o  ^ 

in  (6)  we  get 


a     =  a 
o 


1  =  0^     H 


a 
V.2   "- 


v(v-l)-i(j?+l) 


.-  -.  ^  .^^  =  0,  ;,  0 


This  equation  can  be  satisfied  if  v  =  /+5,    and  from  this, 


+8p  p- 
-     o 


+8p^p- 


(2M)  [(/+5)(i?+2)-i(i+l)]        2{2Ul){2l^+3) 


Therefore 


(7) 


+ 
uj  (X)   = 


+8p  p- 
-     o 


.J?+5 


■^2   x->o  2(2J?+l)(2/+5) 
From  (h)   and  (7)  we  see  that 


J+1 


+  +       +  + 

u:(x)  =  u:  (x)  +  c:  u;  (x)  =    a  ,x^"^ , 

^  ^1     ^  ^2   x->o  ^ 


as  it  can  also  be  inferred  from  the  expansion  of  U-(x)  in  (U.3). 
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